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Abstract 



| We analyze an interacting particle system with a Markov evolution of 

£1+ birth-and-death type. We have shown that a local competition mechanism 

(realized via a density dependent mortality) leads to a globally regular 
behavior of the population in course of the stochastic evolution. 
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ir~i ; 1 Introduction 

Complex systems theory is a quickly growing interdisciplinary area with a very broad 
spectrum of motivations and applications. Having in mind biological applications, 
£f} . S. Levin (see [26] 1 characterized complex adaptive systems by such properties as di- 

versity and individuality of components, localized interactions among components, and 
the outcomes of interactions used for replication or enhancement of components. In 
the study of these systems, proper language and techniques are delivered by the inter- 
acting particle models which form a rich and powerful direction in modern stochastic 
and infinite dimensional analysis. Interacting particle systems have a wide use as 
models in condensed matter physics, chemical kinetics, population biology, ecology 
(individual based models), sociology and economics (agent based models). 

In this paper we consider an individual based model (IBM) in spatial ecology 
introduced by Bolker and Pacala [U [5], Dieckmann and Law [B] (BDLP model). A 
population in this model is represented by a configuration of motionless organisms 
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(plants) located in an infinite habit (an Euclidean space in our considerations). The 
habit is considered to be a continuous space as opposed to discrete spatial lattices 
used in the most of mathematical models of interacting particle systems. We need 
the infinite habit to avoid boundary effects in the population evolution and the latter 
moment is quite similar to the necessity to work in the thermodynamic limit for models 
of statistical physics. Let us also mention a recent paper [2] in which a modification 
of the BDLP model for the case of moving organisms (e.g., branching diffusion of the 
plankton) was considered. 

A general IBM in the plant ecology is a stochastic Markov process in the configu- 
ration space with events comprising birth and death of the configuration points, i.e., 
we are dealing with a birth-and-death process in the continuum. In the particular case 
of the BDLP model, each plant produces seeds independently of others and then these 
seeds are distributed in the space accordingly to a dispersion kernel a + . This part of 
the process may be considered as a kind of the spatial branching. In the same time, 
the model includes also a mortality mechanism. The mortality intensity consists of 
two parts. The first one corresponds to a constant intrinsic mortality value m > s.t. 
any plant dies independently of others after a random time (exponentially distributed 
with parameter m). The second part in the mortality rate is density dependent. The 
latter is expressed in terms of a competition kernel a~ which describes an additional 
mortality rate for any given point of the configuration coming from the rest of the 
population, see Section 3 for the precise description of the model, in particular, (3.6). 
The latter formula gives the heuristic form of the Markov generator in the BDLP 
model. 

Assuming the existence of the corresponding Markov process, we derive in Section 
5 an evolution equation for correlation functions k[ n \ n > 1, of the considered model. 
In HI [S] , [B] this system was called the system of spatial moment equations for plant 
competition and, actually, this system itself was taking as a definition of the dynamics 
in the BDLP model. The mathematical structure of the correlation functions evolution 
equation is close to other well-known hierarchical systems in mathematical physics, 
e.g., BBGKY hierarchy for the Hamiltonian dynamics (see, e.g. !3 ) or the diffusion 
hierarchy for the gradient stochastic dynamic in the continuum (see e.g. 21 ). As in 
all hierarchical chains of equations, we can not expect the explicit form of the solution, 
and even more, the existence problem for these equations is a highly delicate question. 

There is an approximative approach to produce an information about the behavior 
of the solutions to the hierarchical chains. This approach is called the closure procedure 
and consists of the following steps. The first step is to cut all correlation functions of 
the higher orders and the second one is to subscribe the rest correlation functions by 
the properly factorized correlation functions of the lower orders. As result, one obtains 
a finite system of non-linear equations instead of the original linear but infinite system 
of a hierarchical type. This closure procedure is essentially non-unique, see [7]. 

The aim of this paper is to study the moment equations for the BDLP model by 
methods of functional analysis and analysis on the configuration spaces developed in 
|13j . |14j . |15j and already applied to the non-equilibrium birth-and-death type continu- 
ous space stochastic dynamics in [TB], [18]. We obtain some rigorous results concerning 
the existence and properties of the solution for different classes of initial conditions. 
One of the main question we clarify in the paper concerns the role of the competition 
mechanism in the regulation of the spatial structure of an evolving population. More 
precisely, considering the model without competition, i.e., the case a~ = 0, we arrive 
in the situation of the so-called continuous contact model [9], [H], |22] , In the ecolog- 
ical framework, this model describes free growth of a plant population with the given 
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constant mortality. We note that (independently on the value of the mortality m > 0) 
the considered contact model exhibits very strong clustering that is reflected in the 
bound (3.5) on the correlation functions at any moment of time t > 0. Note that this 
effect on the level of the computer simulation was discovered already in [2] and now it 
has the rigorous mathematical formulation and clarification. A direct consequence of 
the competition in the model is the suppression of such clustering. Namely, assuming 
the strong enough competition and the big intrinsic mortality m, we prove the sub- 
Poissonian bound for the solution to the moment equations provided such bound was 
true for the initial state. Moreover, we clarify specific influences of the constant and 
the density dependent mortality intensities separately. More precisely, the big enough 
intrinsic mortality m gives a uniform in time bound for each correlation function and 
the strong competition results ensure the regular spatial distribution of the typical 
configuration for any moment of time that is reflected in the sub-Poissonian bound. 
Joint influence of the intrinsic mortality and the competition leads to the existence 
of the unique invariant measure for our model which is just Dirac measure concen- 
trated on the empty configuration. The latter means that the corresponding stochastic 
evolution of the population is asymptotically exhausting. 

We would like to mention also the work [I(Tj in which the BDLP model was studied 
in the case of the bounded habit in the stochastic analysis framework. The latter case 
differs essentially from the model we consider in the present paper as well as main 
problems studied in |10] , which are related to the scaling limits for the considered 
processes. 

2 General facts and notations 

Let B(R d ) be the family of all Borel sets in R d . B b (R d ) denotes the system of all 
bounded sets in B(R d ). 

The space of n-point configuration is 

To"' = r o?^ : = { V C R d j W = n) , ne No := N U {0}, 

where |^4| denotes the cardinality of the set A. The space := Tq '2 for A € B b (R d ) is 
defined analogously to the space Tq . As a set, Pq™' is equivalent to the symmetrization 
of 

= { (xi, . . . , x n ) G (RT| x k ± Xl if k + l} , 

i.e. (WL d ) n /S n , where S„ is the permutation group over {1, . . . , n}. Hence one can in- 
troduce the corresponding topology and Borel a- algebra, which we denote by 0(Tq ) 
and B(Tq), respectively. Also one can define a measure as an image of the 

product of Lebesgue measures dm(x) = dx on (R d ,S(R d )). 
The space of finite configurations 

To := U P<"' 

nS« 

is equipped with the topology which has structure of disjoint union. Therefore, one 
can define the corresponding Borel cr-algebra B(Tq). 
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A set B 6 B(T Q ) is called bounded if there exists A G B b (R d ) and JV € N such that 
B C |_C=o r A U) • The Lebesgue— Poisson measure A z on To is defined as 

oo n 

A z := > —m y '. 

n=0 

Here z > is the so called activity parameter. The restriction of A z to Fa will be also 
denoted by A z . 

The configuration space 

r := I 7 C R d | |-7 n A| < co, for all A G S b (R d )| 

is equipped with the vague topology. It is a Polish space (see e.g. |15]). The cor- 
responding Borel cr-algebra B(T) is defined as the smallest a-algebra for which all 
mappings Na : T — > No, Na^) := |7 l~l A| are measurable, i.e., 

B(r) =a(iV A |AGS f) (R d )) . 

One can also show that P is the projective limit of the spaces {rA}AGBi,(R d ) w.r.t. the 
projections pa : T —> Ta, PA (7) := 7a, A G B b (R d ). 

The Poisson measure n z on (r, B(P)) is given as the projective limit of the family 
of measures {nt}AeB b (R d )> where is the measure on Pa defined by := e _zm ' A ^A z . 

We will use the following classes of functions: Lp s (Po) is the set of all measurable 
functions on Po which have a local support, i.e. G G //^(Po) if there exists A G S(,(R d ) 
such that G tr \r A = 0; -Bbs(ro) is the set of bounded measurable functions with 
bounded support, i.e. G \t \b~ for some bounded B G B(To)- 

On r we consider the set of cylinder functions TL° (P) , i.e. the set of all measurable 
functions G on (P,S(r))) which are measurable w.r.t. Sa(P) for some A G Sf,(R ). 
These functions are characterized by the following relation: ^(7) = F fr A (ta)- 

The following mapping between functions on To, e.g. L? s (Po), and functions on V, 
e.g. TL°(T), plays the key role in our further considerations: 

where G G L^ro), see e.g. [131 1241 [25] , The summation in the latter expression is 
taken over all finite subconfigurations of 7, which is denoted by the symbol 77 <g 7. 
The mapping K is linear, positivity preserving, and invertible, with 

K-'F^^^-V^HO, ver t) . (2.2) 

Sen 

Let .A4f m (r) be the set of all probability measures /i on (P, S(P)) which have finite 
local moments of all orders, i.e. J r ta |" A*(c?7-) < +00 for all A G Sb(R d ) and n G No- 
A measure p on (r , B(r )) is called locally finite iff p(A) < 00 for all bounded sets A 
from B(To). The set of such measures is denoted by A^if(Po). 

One can define a transform K* : _Mf m (P) — » A^if(ro), which is dual to the K- 
transform, i.e., for every fi G Mf m (T), G G Bbs(Po) we have 



f KG(f)n(<h,) = [ G{ n )(K*n){dri). 
Jr Jr 
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The measure p^, := K* p is called the correlation measure of fi. 

As shown in |13] for p £ .Mf m (r) and any G £ L 1 (ro,p, 1 ) the series (|2.ip is p-a.s. 
absolutely convergent. Furthermore, KG £ L 1 (F,/i) and 

G{ V ) pMv) = f(KG)(j) ii(dry). (2.3) 

A measure p £ A^f m (r) is called locally absolutely continuous w.r.t. n z iff pa '■= 
pop^ 1 is absolutely continuous with respect to 7r^ for all A £ 6a (R d ). In this case 
pp. := K* p is absolutely continuous w.r.t \ z . We denote 



The functions 



-„(«) 



(2.4) 



i.(«)( Tl T \ ._ *v({a?i,... ,£„}), if (an,. ..,£„) £ (R d )" 
0, otherwise 



are the correlation functions well known in statistical physics, see e.g [28], [29] , 

We recall now the so-called Minlos lemma which plays very important role in our 
calculations (cf., |20j). 

Lemma 2.1. Let n £ N, n > 2, and z > be given. Then 

/ .../ G(7?i U .. . U r) n )H(rn, . . . ,ri n )d\z(rn) .. -dX^^n) 



L 



G(n) J2 H{ Vll ...,r, n )d\M 



for all measurable functions G : To i— > R and H : To X . . . X To i— > R mift respect to 
which both sides of the equality make sense. Here Vnij]) denotes the set of all ordered 
partitions of n in n parts, which may be empty. 



3 Description of the model 

In the present paper we study the special case of the general birth-and-death pro- 
cesses in continuum. The spatial birth-and-death processes describe evolution of con- 
figurations in R d , in which points of configurations (particles, individuals, elements) 
randomly appear (born) and disappear (die) in the space. Among all birth-and-death 
processes we will distinguish those in which new particles appear only from existing 
ones. These processes correspond to the models of the spatial ecology. 

The simplest example of such processes is the so-called "free growth" dynamics. 
During this stochastic evolution the points of configuration independently create new 
ones distributed in the space according to a dispersion probability kernel < a + £ 
L 1 (R d ) which is an even function. Any existing point has an infinite life time, i.e. 
they do not die. Heuristically, the Markov pre-generator of this birth process has the 
following form: 

(L+F)( 7 ) = x + ]T / a + (x - y)D+F(rr)dx, 



6 



D. Finkelshtein, Yu. Kondratiev, O. Kutoviy 



where 

DiF(j) = F(jUx)-F(j), 

and >c + > is some positive constant. 

The existence of the process associated with L+ can be shown using the same 
technique as in [9], [22]. Let fit be the corresponding evolution of measures in time 
on .A4f m (r). By k\ , n > we denote the dynamics of the corresponding n-th order 
correlation functions (provided they exist). Note, that each of such functions describes 
the density of the system at the moment t. 

Then, using (|2.3|l . for any continuous <p on R d with bounded support, we obtain 

difd V^)^ ( x > dx = ^ J (Vj7>dA**(7) = J L +^,~l)dM"i) 

Jr Jm d 
= x + / (y3(s)(a + * k$P)(x)dx, 

where * denotes the classical convolution on R d . Hence, k\ grows exponentially in 
t. In particular, for the translation invariant case one has k { 1] (x) = 4 1 ' > and as a 
result 

fc'^e^C- (3-D 

One of the possibilities to prevent the density growth of the system is to include 
the death mechanism. The simplest one is described by the independent death rate 
(mortality) m > 0. This means that any element of a population has an independent 
exponentially distributed with parameter m random life time. The independent death 
together with the independent creation of new particles by already existing ones de- 
scribe the so-called contact model in the continuum, see e.g. |22| . The pre-generator 
of such model is given by the following expression: 

(£cmF)( 7 ) = m^fl.-F(7) + (L+F)( 7 ) 

= mY J DZF( n ) + x + Y J ( a+(x-y)DtF(>y)dx, 

where 

D-F( 1 ) = F( 1 \x)-F( 7 ). 

The Markov process associated with the generator Lcm was constructed in [22] , 
This construction was generalized in [9] for more general classes of functions a + . Let 
us note, that the contact model in the continuum may be used in the epidemiology 
to model the infection spreading process. The values of this process represent the 
states of the infected population. This is analog of the contact process on a lattice. Of 
course, such interpretation is not in the spatial ecology concept. On the other hand, 
contact process is a spatial branching process with a given mortality rate. 

The dynamics of correlation functions in the contact model was considered in [17] . 
Namely, taking m = 1 for correctness, we have for any n > 1, t > the correlation 



IBM Model with Competition in Spatial Ecology 



7 



function of n-th order has the following form 



k t [xi, . . . , Xn) 



(^1 , • ■ • , X n ) 



(3.2) 



/ 6 

■A) 



)(*-«) 



('-»)i' 



where 



- y) |V n) (a;i, ...,Xi- 



i,y,xi+i, 



x n ) — k {n \xi, . . . ,x„) \ dy 



and the symbol Xi means that the i-th coordinate is omitted. Note that L z a+ is a 
Markov generator and the corresponding semigroup (in L°° space) preserves positivity. 
It was also shown in [17], that if there exists a constant C > (independent of n) such 
that for any n > and (xi, . . . , x„) G (R d ) n 

ko(xi, . . . ,x n ) < n\ C™, 

then for any t > and almost all (a. a.) (x\, . . . , x n ) G (R d ) n (w.r.t. Lebesgue measure) 
the following estimate holds for all n > 



k { t n) ( Xl , ...,x n )< /(t)"(l + a ) n e n( " + - 1)t (C + t) n n\ 



(3-3) 



Here 



a = a 



x + (t) := max ^1, x + , x + e 



-l)t 



For the translation invariant case the value x — 1 is critical. Namely, from (|3.2|l 
we deduce that 

fcW = e^- 1 *^. (3.4) 

Therefore, for x + < 1 the density will exponentially decrease to (as t — * oo), for 
k + > 1 the density will exponentially increase to oo, and for x + = 1 the density will 
be a constant. One can easily see from the estimate (13. 3p that, in the case k + < 1, 
the correlation functions of all orders decrease to as t — > oo. On the other hand, 
for fixed t, the estimate (|3.3|l implies factorial bound in n for kf . As result, we may 
expect the clustering of our system. To show clustering we start from the Poisson 
distribution of particles and obtain an estimate from below for the time evolutions of 
correlations between particles in a small region. 

Hence, let x+ < 1, k { n) = C n . Let B is some bounded domain of R d such that 

a := inf a + (x — y) > 0. 

x,y£B 

Let P — min{ax + ,C}. For any {x^,X2} C B, formula (|3.2[) implies 

l ( 2 ) I \ \ + f 2(x+-l)(t-s) j \ r,£)2, 2(xr+-l)t 

fe t (xi,X2) > 2Cx a e A 'as > 2p te K ' . 
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We consider t > 1. One can prove by induction that for any {x\, . . . , x n } C B, n > 2 
^(ii,...,^)^^"- 1 ^! (3.5) 



Indeed, for n — 2 this statement has been proved. Suppose that (|3.5[) holds for n — 1. 
Then, by (|3.2p . one has 



fc t (n) (x!, ...,x n )>x + e n( " - 1)(t - s) n/3 (n - 1) e (n - 1)(,c 1)s (n - l)!( n - l)ada 
Jo 

\ £>n i n(x~^~ — l)t I — — l)s t ^ r>n n(x~^~ — l)t i 

> p n'.e v ' e ds > p e ' n\ . 

Jo 

As it was mentioned before, the later bound shows the clustering in the contact model. 
All previous consideration may be extended for the case m / 1: we should only replace 
1 by m in the previous calculations. 

As a conclusion we have: the presence of mortality (m > x + ) in the free growth 
model prevents the growth of density, i. e. the correlation functions of all orders decay 
in time. But it doesn't influence on the clustering in the system. One of the possibilities 
to prevent such clustering is to consider the so-called density dependent death rate. 
Namely, let us consider the following pre-generator: 



m + yt a (x — y) 



(£F)(7) = E 

Here < a~ £ L 1 (R d ) is an arbitrary, even function such that 

a~ (x)dx = 1 



D~F(i) (3.6) 



/ 



(in other words, a is a probability density) and >t > is some positive constant. It 
is easy to see that the operator L is well-defined, for example, on TL°(T). 

The generator (|3.6p describes the Bolker — Dieckmann — Law — Pacala (BDLP) model, 
which was introduced in [UEIE]. During the corresponding stochastic evolution the 
birth of individuals occurs independently and the death is ruled not only by the global 
regulation (mortality) but also by the local regulation with the kernel x~ a~ . This 
regulation may be described as a competition (e.g., for resources) between individuals 
in the population. 

The main result of this article is presented in Section [5] Theorem 15.11 It may be 
informally stated in the following way: 

If the mortality m and the competition kernel x~ a~ are large enough, then the 
dynamics of correlation functions associated with the pre-generator ()3.6|) preserves 
(sub-)Poissonian bound for correlation functions for all times. 
In particular, it prevents clustering in the model. 

In the next sections we explain how to prove this fact. In the last section of the 
present paper we discuss the necessity to consider "large enough" death. 
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4 Semigroup for the symbol of the generator 

The problem of the construction of the corresponding process in F concerns the pos- 
sibility to construct the semigroup associated with L. This semigroup determines the 
solution to the Kolmogorov equation, which formally (only in the sense of action of 
operator) has the following form: 

lit= LF *> F Lo= F °- (4- 1 ) 

To show that L is a generator of a semigroup in some reasonable functional spaces 
on F seems to be a difficult problem. This difficulty is hidden in the complex structure 
of the non-linear infinite dimensional space F. 

In various applications the evolution of the corresponding correlation functions (or 
measures) helps already to understand the behavior of the process and gives candidates 
for invariant states. The evolution of correlation functions of the process is related 
heuristically to the evolution of states of our IPS. The latter evolution is formally given 
as a solution to the dual Kolmogorov equation (Fokker — Planck equation): 

— = L pit, /it \ t=Q = fio, (4.2) 

where L* is the adjoint operator to L on A4} ln (F), provided, of course, that it exists. 

In the recent paper [16], the authors proposed the analytic approach for the con- 
struction of a non-equilibrium process on F, which uses deeply the harmonic analysis 
technique. In the present paper we follow the scheme proposed in [16] in order to con- 
struct the evolution of correlation functions. The existence problem for the evolution 
of states in Mj ul (F) and, as a result, of the corresponding process on F is not realized 
in this paper. It seems to be a very technical question and remains open. 

Following the general scheme, first we should construct the evolution of functions 
which corresponds to the symbol (K-image) L — K~ 1 LK of the operator L in L 1 - 
space on To w.r.t. the weighted Lebesgue — Poisson measure. This weight is crucial 
for the corresponding evolution of correlation functions. It determines the growth of 
correlation functions in time and space. Below we start the detailed realization of the 
discussed scheme. 

Let us set for n 6 To 

where a* denotes either a~ or a + . 

Proposition 4.1. The image of L under the K-transform (or symbol of the operator 
L) on functions G £ -Bi, s (Fo) has the following form 

(LG)( V ) := (K~ x LKG) (77) 

m\n\ + x'E a {n)^j G{n) - x~ Yl a ~ i x ~ v) G ( n \v) 



+ x+ I J2 a+ ( x ~ y) G ^ \ y) u x ) dx 

Jmd yen 

+ x + I y^ j a + {x- y)G{r) U x)dx. 
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For the proof see [8]. 

With the help of Proposition 14.11 we derive the evolution equation for quasi- 
observables (functions on To) corresponding to the Kolmogorov equation (|4.ip . It 
has the following form 

^=LG t , Gt 1^= Go. (4.3) 

Then in the way analogous to those in which the corresponding Fokker-Planck equa- 
tion (|4.2[) was determined for (|4.1[) we get the evolution equation for the correlation 
functions corresponding to the equation (|4.3[1 : 

dkt ~ 

— = L*k t , h\ t=0 =k . (4.4) 

The precise form of the adjoint operator L* will be given in Section [5] It is very 
important to emphasize that in the papers [4j [5] the equation (|4.4[) was obtained from 
quit heuristic arguments and, moreover, it was considered as the definition for the 
evolution of the BDLP model. 

Let A be the Lebesgue-Poisson measure on To with activity parameter equal to 1. 

For arbitrary and fixed C > we consider the operator L as a pre-generator of a 
semigroup in the functional space 

C c :=L\Y ,C M \{dn)). (4.5) 

In this section, symbol stands for the norm of the space (|4.5[) . 

For any u> > we introduce the set Ti(u>, 0) of all densely defined closed operators 
T on Cc, the resolvent set p(T) of which contains the sector 

Sect :={£eC | | arg C| < | + u/f , to > 

and for any e > 

IKT-ar 1 !! < |argC|<!+^- £ , 

where M E does not depend on £. 

Let TL{ijj, 6), 6 € K denotes the set of all operators of the form T = To + 6 with 
T,)G%0). 

Remark 4.1. It is well-known (see e.g., that any T £ 1~L{ui, 9) is a generator of 

a semigroup U(t) which is holomorphic in the sector | arg t\ < uj. The function U(t) 
is not necessarily uniformly bounded, but it is quasi-bounded, i.e. 

\\U(t)\\ < const|e 9t | 

in any sector of the form \ arg t\ < u — e. 

Proposition 4.2. For any C > 0, m > 0, and x~ > the operator 
(L G)( V ) := - (mM + (r?)) G(r?), 

D(L ) = {GeCc\ (m\v\ + * E a ~ (77)) G( v ) e C c } 
is a generator of a contraction semigroup on Cc- Moreover, Lo £ Tt(uj,0) for all 

we (o,|). 
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Proof. It is not difficult to show that Lq is a densely defined and closed operator in 
Cc 

Let < lu < § be arbitrary and fixed. Clear, that for all £ € Sect (■§ + 

\m\ri\ + x~E a (77) + C| > 0, r; 6 IV 

Therefore, for any C, £ Sect + a;) the inverse operator (Lo — , the action of 

which is given by 



[{L -ar 1 G]{r,) = - 



m\rj\+>c-E a (77) + C 



G(v), 



(4.6) 



^, ifReC<0, 



(4.7) 



is well defined on the whole space Cc- Moreover, it is a bounded operator in this space 
and 

( jij, ifReC>0, 

\\(L - am < i 

where the constant M does not depend on £. 

The case Re f > is a direct consequence of (|4.6[) and inequality 

m\r]\ + K~E a ~ {n) +RcC > Re C > 0. 

We prove now the bound (|4.7|l in the case Re£ < 0. Using (|4.6[) . we have 

1 



|(Lo -Cl) _1 G||c = 



rn\ ■ I +x-E a (0 + C 
ICI 



G(- 



m| ■ I + x-E a (■) + C, 



Since £ € Sect (f + w), 



Hence, 



|ImC| > |C||sin(^-+w)| = |C|cosw. 



■G(. 



ICI 



< jcl < 1 



+ K~E a (??) + C| |I m CI cosw 



M 



and (53 is fulfilled. 

The rest of the statement of the lemma follows directly from the theorem of Hille — 
Yosida (see e.g., PU). □ 

We define now 

(LiG)(i7):=x-£; J2 a-(x-y)G{r,\y), G £ D(Li) := D(Lo). 
The lemma below implies that the operator L 1 is well-defined. 

Lemma 4.3. For any S > there exists Co := Co(S) > such that for all C < Co 
the following estimate holds 



\LiG\\c<a\\L G\\ c , G € D(Li) 



(4.8) 



with a = a(C) < S. 
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Proof. By modulus property 

||LiG||o= / |LiG(77)|C W A(dr;) (4.9) 

can be estimated by 

x~ f E a~(x-y)\G(v\y)\C M \(d V ). (4.10) 
-' r o xev yev\x 

By Minlos lemma, (|4.10|l is equal to 



Jr Jr* x€v 



x \ C\n\\G(j])\C w \(dn) < C\\L G\\ C . 

r ~ m 



Therefore, (|4.8|) holds with 

x~C 

a — . 

m 

Clear, that taking 

G = ^ 
x 

we obtain that a < S for C < Co- □ 
We set now 

{L 2 G){n) := {L 2 ^ + G){n) = x + / ^a+tx-rfGf^rfUx)^ 

jRd yen 

GED(L 2 ) := D(L ). 

The operator (L 2 , D (L2)) is well defined due to the lemma below: 

Lemma 4.4. For any 8 > there exists x^ := Xg(8) > such that for all x + < Xq 
the following estimate holds 

\\L 2 G\\ C <a\\L G\\ c , G G D(L 2 ), (4.11) 

with a — a(x + ) < 5. 

Proof. Analogously to the previous lemma we estimate 

\\L 2 G\\o= [ \L 2 G(n)\C M X(dr,) (4.12) 



by 

x + [ [ J2 a+ ( x -y)\G{(ri\y)Ux)\C M dx\(dn). (4.13) 
Jr Jm* yev 

By Minlos lemma, (|4.13[) is equal to 

^ I E / a + (x-y)\G( V )\C^dyX(dr,)<— \\L G\\c. 
Taking Xq = 8m we prove the lemma. □ 
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The operator defined as: 

(NG)(r,) = \r,\G(r,), G 6 D(L ) (4.14) 
is called the number operator. 

Remark 4.2. We proved, in particular, that for G G D(Lo) = D(Li) = D(L,2) 

\\LiG\\ c < x-C\\NG\\ c , 
\\L 2 G\\ C < x + \\NG\\ c . 

Finally, we consider the last part of the operator L: 

{L 3 G)( V ) :=x + [ ^fl+fi-^U^i, D(L 3 ) := D(L ). 

Jmd yev 

Lemma 4.5. For any 8 > and any x + > 0, C > such that 

x + E a+ (rj) < SC (n~E a ~ (rj) + m\r)\) (4.15) 
the following estimate holds 

||i3G||o<o||L G||c, G G D(L 3 ), (4.16) 

with a = a(x + , C) < 5. 

Proof. Using the same tricks as in the two previous lemmas we have 
\\L 3 G\\ C = [ |L 3 G(7 ? )|C l ' j| A(dr ? ) 



< h + f [ ^ j a + {x-y)\G{r)yjx)\C M dx\{dri). (4.17) 
Jt V ySt) 

By Minlos lemma, (|4.17p is equal to 

Er + {r,)\G(r,)\cW\{d V ). 

The assertion of the lemma is now trivial. □ 

Theorem 4.6. Assume that the functions a~ , a + and the constants x~ , x + > 0, 
m > and C > satisfy 

Cx~aT > 2x + a + , (4.18) 
m > 2 (x"C + x + ) . 
Then, the operator L is a generator of a holomorphic semigroup U t , t > in Cc- 



Proof. The statement of the theorem follows directly from Remark 14.21 Lemma 14.51 
and the theorem about the perturbation of holomorphic semigroup (see, e.g. [12]). 
For the reader's convenience, below we give its formulation: 

For any T £ TI(uj, 0) and for any e > there exists positive constants a, 8 such 
that if the operator A satisfies 

\\Au\\ < a\\Tu\\+b\\u\\, u e D(T) C D(A), 
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with a < 8, b < 5, then T + A <= H(u - e, a). 

In particular, if 9 = and 6 = 0, then T + A G H,(lu — e, 0) 

Following the proof of this theorem (see, e.g. [12]) and taking into account the 
fact that Lo G H(lo, 0) for any w 6 (0, |), one can conclude in our case that 8 can be 
chosen equal to | . This is exactly the reason of appearing multiplicand 2 at the l.h.s. 
of (|4T8j) . " □ 



5 Evolution of correlation functions 

Let us consider the evolution equation (|4.4|l . which corresponds to the operator L* 

^j-=L*k t , h\ t=0 =k . 
Using the general scheme, proposed in [S] we find the precise form of L*: 

L* k(rj) — — ^m|r7| + x E a (rf)^ k(-n) + x + a + (x — y)k{rj \ x) 

xEr} y£r)\x 

+ J2 a+ ( x ~y) k (( r i\y)ux)dx 

-*-[ r«- (l -# ( ,u#. 

The main questions which we would like to study now are the existence and properties 
of the solution to the hierarchical system of equations (|4.4|l . The answers to these 
questions are given in the following theorem 

Theorem 5.1. Suppose that all assumptions of Theorem \4- 6] are fulfilled. Then for 
any initial function ko from the class 

Kc ■= |fc:r ^R I k-C~ M GL°°(r ,A)j 

the corresponding solution kt to exists and will be again the function from K-c 

for any moment of time t > 0. 

Proof. Following the scheme proposed in |16| . we construct the corresponding 
evolution of the locally finite measures on To. In order to realize this construction 
we consider the dual space Kc to the Banach space Cc- The duality is given by the 
following expression 

{{G,k)):= I G-kd\, GeCc- (5.1) 
J v 

It is clear that Kc is the Banach space with the norm 

||fc|| := ||C-l-lfc(-)|U-(r ,A). 
Note also, that k ■ C~''' G L°°(ro, A) means that the function k satisfies the bound 
\k{rj)\ < constc' 11 ' for A-a.a. r\ G Tq. 
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The evolution on JCc, which corresponds to Ut, t>o constructed in Theorem 14.61 may 
be determined in the following way: 

«G, fct» := ({UtG, fco)) • 

We denote 

Utk ■— k t . 

Using the same arguments as in [16], it becomes clear that kt = ?7 t *fco is the solution 
to (|4.4[1 in the Banach space K-c- □ 
It is important to emphasize that in the case of a~ = and x + < m 

kl —* 0, t — y 0, for any n > 1, 

see e.g. |17j . Therefore, we may expect that the correlation functions of our model 
satisfy this property as well. 

6 Stationary equation for the system of correla- 
tion functions 

Let us consider for any a £ K the following Banach subspace of KLc'- 

/Cg := {kelCc |fc (o) (0) = a}. 

In this section we study the existence problem for the solutions to the stationary 
equation 

L*fc = (6.1) 
in K-c- The main result is formulated in the following way: 



Theorem 6.1. Suppose that 



^ + ^<1 (6.2) 
m m C 



x a > x a 

then the solution k = (fe'™^) n >o to 16. lp is unique in K} c and such that 

fc (n) =0, n > 1. 

Proof. Let 

(L*k) (77) = 0. 

The latter means that 

( m \n\ + x~ E a (77) J k (77) = — x~ ~S~] / k (y U 77) a~ (x — y) dy + 

+x+ J2 a + {x~ y)k{n \ x) + x+ J a+ ( x ~ v) k ((.V \ ^ u s ) d:r - 

The last relation holds for any k G /C^ at the point 77 = 0. Hence, one can consider it 
on K° c . 
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Let us denote for 77 7^ 

(Sk)(r,) = * V / fe (y U 17) a" (a; - y) + 



+ 



m 


\v\ 


+ x- 


E a ~ 


(v) 












m 


\v\ 






(v) 












m 


[>l\ 


+ x- 


E a ~ 


(v) 



Yl Yl a+ i x -y) k (v\ x ) + 

x^rj y£r}\x 

f 5>+(,-y M ,\y)u^ 



and 
Let 

then 



(Sk) (0) = 0. 

!|fc|| c =esssupM : 
ver oiii 



l|Sfc|| c 

< ll&llri ess sup — ■ — : — — — — — — — V f a~(x — y)dy 

%er o \{0} m\n\ + x~E a (77) 

+ ± -FT- esssup — — _ - > > a + (a;-y) 

+ ||fc|| c esssup / \^a + (x — y)dx 

ver \{9} m \rj\ + E a (rj) J M d ^ 



,,,,, x ,,,,, 1 „,,, fCx x + 1 

— + Wc- + WI ^ = W c — + 



if 

As result, 



x + £ a+ (77) < x E a {ji) + m\<q\ 



\\s\\ < -c>c- + ->c + + 1. 

m m C 



The assertion of the theorem is now obvious. □ 



Remark 6.1. For any C > 1 one may chose x~ > and m > suc/i t/iar; Iff.gp is 
satisfied. The latter means, that, asymptotically, our system exhausted to the system 
with the stationary state S^^y) (the Dirac measure concentrated on the empty con- 
figuration In other words, the population evolving due to the BDLP dynamics is 
asymptotically degenerated. 

7 Further developments 

In Theorem 15.11 we have shown that functions k t is bounded by C™ for all t > 0, 
provided that ko satisfies initially the bound of the same type. Using approximation 
arguments (see e.g. [16], [18]) one may prove that the corresponding time evolution of 
the correlation function will be also correlation function for some probability measure 
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on r. We suppose to discuss this problem as well as other probabilistic aspects of 
the BDLP model in a forthcoming paper. The main aim of the present paper is to 
analyze evolution of correlation functions. Namely, we have shown that dynamics of 
correlation functions stays in the space ICc- This property seems to be very strong. To 
show that system of correlation functions evolving in time stays in the same space is 
already difficult even for the contact model. Namely, f|3 . 3|) implies that the evolution 
of correlation functions at some moment of time t may leave the space 



The reason is that C may depend on t, which is true at least for the case x > 1 
(m = 1 at the moment). Hence, we may expect that the dynamics of correlation 
functions for the contact process lives in some bigger space. Of course, this is possible 
only for x + < 1 since for x + > 1 density tends to infinity. Hence, let us consider the 
case >c + = 1. One candidate for such bigger space is 



Note, that the invariant measure of the contact process belongs to this space (see |17l 
Theorem 4.2]), provided that d > 3, a + has finite second moment w.r.t. the Lebesgue 
measure and the Fourier transform of a + is integrable on R d . Below we show that the 
evolution of correlation functions at any moment of time t is a function from IZc- 

Indeed, let x + — 1 and suppose that there exists C > such that for any n > 1 
and for any x\ , . . . , x n G R d 



Then, it is clear that k G Tie- Now, suppose that k { " 15 < C" -1 ((n- 1)!) 2 . We prove 
the corresponding inequality for fcj using the mathematical induction. By (|3.3[) we 
have for any xi, ■ . . , x n £ R d 





k {xi, . . . , 



xn) < ic"(n!) 2 . 



(7.1) 



+ 




n n 



Xi — Xj)da 




where for / G L 1 (R d ) 




For the bound above we have used the fact, that for any 1 < i ^ j < n 



Ll + a + (x t -Xj) =L J a 



a + (xi - Xj) = (L a +a + )(xi - Xj), Xi,Xj G R d . 



This relation can be easily checked by simple computations. 
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Note, that L a + is a generator of the Markov semigroup which preserves positivity 
in L 1 (R d ). Hence, 

g t {x):= [ (e 2(t ~ s)L -+a + ^j (x)ds>0, x£~R d ,t>0, 



and g t G L (Mr). Then we have 

9t{x) = \gt(x)\ = — / e px g t (p)dp 
(2n) d J Rd 



\ Z7T ) JR d JO 



where symbol / denotes the Fourier transform of the function / G L (R ). Therefore, 
1 f 1 - e 2t(a+(p)_1) i If \a + (p)\ 

9t{x) ~ Wr L 2 ( i-a + (p)) '° (p)|rfp * w L ir^j*- 

It was shown in 17 that under the conditions posed on function a + for the case of 
invariant measure 

D:=[ -K^d P <oc. 
J R d 1 - a+(p) 

Finally, if additionally 

C> D 



then we obtain from (|7.ip 

kl n) ( Xl , ...,*»)< ic"(n!) 2 + ic"- 1 ^ - l)!) 2 n(n - 1)^ < G"(n!) 2 . 

As result, kt G 7?.c for all t > 0. 

Therefore, the dynamics of correlation functions for the contact model stays in 
hence, this dynamics is really very clustering for >c + = m = 1. As before, we may 
extend our consideration on the case ra/1. 

Summarizing previous results in this section we claim that the presence of the big 
mortality and the big competition kernel prevents clustering in the system making 
it sub-Poissonian distributed. But, is it really necessary to add "big" mortality and 
competition kernel? Below we discuss this problem. 

If we want to study the quasibounded semigroup with the generator L on Cc for 
some C > then, naturally, this generator should be an accretive operator in Cc- 
Hence, for some b > the following bound should be true 

jf sgn (G (r,)) ■ (jti — 61)G) (rj) d\ c (ij) < 0, V G G D(L), 

since 

C^^dX(rj) — d\c{rf). 
Let us define the "diagonal" part of the operator L: 



{L dlag G) (t?) := -m|ij|G(ij) - k~ E a ( V )G( V ) + x+ f V\ 

Ju " yen 

c n > 1 

G= (o,0,G (n) ,0,o) , G (n) G L 1 ((R d ) n ). 



;a + (s-y)G((7i\2/)Ux)da; 

and consider for some n > 1 
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Then 



x + I E a + (x-y)G {n) (r)Ux)dx, \r)\ = n - 1 
-xT E E «-(i-s)G W (t)\!/), |»?|=n + l 

0, otherwise 



Note that sgn (G (77)) = if |?7| / n. 
Therefore, for arbitrary n > 1 



> /„ := 

'r 



jT sgn (G (r?)) • ((L - 61) G) (77) d\ c (tj) 
= /" sgn (G (r?)) ■ [(L diag - 61) G M ) (77) dA c (r?) 
= ^ ^ sgn(G (n) (*<">)) ((L dm9 - 6l)G (n) ) (x (n) )dx ( 



Let us fix some t > and A £ /3i,(R d ). Set for n > 1 

n 

G W (x W ) = T ]J XA (a*) = ri rW ({x ( ">}) € L 1 ((R d ) n ). 
k=l A 

Then, the equality 

n 

sgn(G (n) ( a ; ( " ) )) = nXA(x fc ) 
fc=i 

implies 

= / -mn J] XA (a:*) - x"£ a " (x M ) XA (a*) 

V fc=i fc=i 

V* a + (x - a;^) IT xa (ifc) xa (a;) dx dx (n) - 6 / IT xa (xk) dx 
= — x~ I E a (x' n ')dx' n ' + x + TT / dxfc / / a + (x — Xj)dxdxj 

J A™ j=1 J A"- 1 JaJA 

- (6 + mn) \A\ n 

= -*:-[ £ a ~(x (n) )dx (n) +x + n|A|™- 1 / / a+(x - y)dxdy - (6 + ttiti) |A| n • 

J A n J A J A 

We suppose, in fact, that for any n > 1 



(n) 



20 



D. Finkelshtein, Yu. Kondratiev, O. Kutoviy 



Since E a (rj) = for |r?| < 1 we get 

°° °° in s^vn °° in/in r 

> E '» = ~m £ nL£- | A|» - x- £ *£- / (*<»>)*,<»> 

n=l n=l ' n=l ' J A ™ 

+* + £^t-»i a i b ~7 / « + (^y)^-bE^ri A r 

= -mtC|A|e Ct|A| -x" / E a ~ (n) d\ C t {n) + x+ Cte CtW [ [a + (x-y)dxdy 
Jr A J a J a 

- b (e ct ' A ' - l) 

= -mtC\A\e cm - x~C 2 t 2 f [ f a~ (x - y) dxdydXct (r?) 

ir A 7a Ja 

+ x+Cte Ct|A| | J a + {x-y)dxdy-b(e ctw -l) 

_ e ct|A| ^ I x + / / a +^ x _ y^dxdy — x~Ct I I a~ (x — y) dxdy — m\A\ 
L V Jaja Jaja 

-b(l-e- ct ^)' . 
Therefore, for any t > and any A £ Si,(R d ) 

0>x+[ [a+(x- y)dxdy - x~ Ct [ [ a~ (x - y) dxdy - m |A| 
7a 7a Jaja 

( 1 _ e ~Ct\A\\ 

- b K ct 

Suppose that there exists z > such that 

a+ i x ) > za ~ ( x ) i * £ R d i 

then taking for some e > 



zx 

t = e—->0 

x C 



we obtain 
B > 

bx 



(1 — e)x + z I a (x — y)dxdy — in |A| 

JAJA 

(l-e*p(-g£|A|)) ~ ((l- E )x + *-m)|A|, AT 



ezx 

which contradicts to B < 0. As result, m can not be arbitrary small. 
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